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ABSTRACT 

In 1997 Soker laid out a framework for understanding the formation and shap- 
ing of planetary nebulae (PN). Starting from the assumption that non-spherical 
PN cannot be formed by single stars, he linked PN morphologies to the binary 
mechanisms that may have formed them, basing these connections almost en- 
tirely on observational arguments. In light of the last decade of discovery in the 
field of PN, we revise this framework, which, although simplistic, can still serve 
as a benchmark against which to test theories of PN origin and shaping. 

Within the framework, we revisit the role of planets in shaping PN. Soker 
invoked a planetary role in shaping PN because there are not enough close bi- 
naries to shape the large fraction of non-spherical PN. In this paper we adopt a 
model whereby only ~20% of all 1 — 8 M stars make a PN. This reduces the 
need for planetary shaping. Through a propagation of percentages argument, 
and starting from the assumption that planets can only shape mildly elliptical 
PN, we conclude, like in Soker, that ~20% of all PN were shaped via planetary 
and other substellar interactions but we add that this corresponds to only ~5% 
of all 1 — 8 M stars. This may be in line with findings of planets around main 
sequence stars. PN shaping by planets is made plausible by the recent discovery 
of planets that have survived interactions with red giant branch (RGB) stars. 

Finally, we conclude that of the ~80% of 1 — 8 M stars that do not make 
a PN, about one quarter do not even ascend the AGB due to interactions with 
stellar and substellar companions, while three quarters ascend the AGB but do 
not make a PN. Once these stars leave the AGB they evolve normally and can 
be confused with post-RGB, extreme horizontal branch stars. We propose tests 
to identify them. 

Subject headings: planetary nebulae: general — stars: AGB and post-AGB - 
stars: mass loss 
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Introduction 



Sokerl ( 119971 ) established an observationally-based framework whereby planetary nebulae 
(PN) with different morphologies were mapped to the shaping mechanism that could produce 
them. These mechanisms ranged from single stars to a multitude of binary interactions with 
stellar and substellar (brown dwarf and planetary) companions. 

Many discoveries, both observational and theoretical have been made in the last ten 
years that have thicken ed the debate of what shapes PN. Key t o the debate have been 
th e theoretical work of Sokerl (l2006bh and iNordhaus et al.l (120071 ) (but see previous work 
by ISoker fc Zoabil (120021 )). who described the difficulty with which single asymptotic giant 
branch (AGB) stars can sustain rotation and global magnetic fields for long enough to affect 
the geometry of the mass-loss and the shape of the subsequent PN. In the case of single 
stars, rotation and global magnetic fields, have been the leadin g mechanism suggested for the 
shaping of non-spherical PN (e.g., iGarcfa-Segura et al.l 120051 ). Lacking a full understanding 



of how single AGB stars can produce winds that greatly diverge from a spherical distribution, 
sev eral authors star ted looking more favourably to binary origin explanations (for a review 



sec 



De Marcd 12009} ) . where the companion is either a AGB wind shaping agent (i.e., 



via 



gravity), or it induces primary envelope rotation and magnetic fields. 

Irrespective of what camp of the debate each researcher favours, the scheme of 
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tional and theoretical results, there is scope for a revision of that framework. In particular, 
the role of planets in shaping PN was already discussed in the late 90s, but could not be 
observationally quantified at that time. Today, we know more on the statistics of planetary 
systems, and, although many questions still remain, this is enough to update our estimate 
of the role of planets in shaping PN. 

Finally, new planetary discoveries are changing our understanding of stars. New ques- 
tions are being asked as to the influence planets have on stellar evolution, and, conversely, 
how stars affect planet evolution and survival (see for instance the proceedings of the con- 
ference Planets Beyond the Mains Sequence; Bamberg, Germany, August 11-13, 2010). The 
role of planets in PN shaping is topical in this context. 



In § [2] we explain the PN shaping framework laid out by ISokerl (119971 ) . In § [3] we discuss 
ne w resu l ts in t he field of PN (§ !3.ip and planets (§ 13.21) . which allow us to reassess the work 
of ISokerl (119971 ) in § |U In § [5] we present some arguments regarding the detection of naked 
central stars, while in § [6] we conclude. 
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2. The PN shaping framework in 1997 

Soker (1997) pre sent ed a classification of 458 PN , based on the morphological systems of 



Schwarz et al.l (119931 ) and lCorradi &: Schwara (119951). with th e aim of distinguishing between 



the role of stellar and substellar companions. ISokerl ( 119971 ) further classified elliptical PN 
into those with large and small departure from sphericity. He divided PN morphologies into 
four categories and flanked each category with the physical process most likely to give rise 
to that shape using a series of observational arguments. 

(1) Spherical PN are formed by progenitors that did not have a companion, or did not 
interact with their companion. Some departure from sphericity (though no axi-symmetry) 
may be expected if the companion separation is large, but still interacting in some way with 
the primary and the PN is formed over a span of time comparable with the orbital period. 

(2) Bipolar PN are formed by a close stellar companion that avoided a common envelope 
phase, or that entered the common envelope only late in the evolution. The stellar companion 
accretes a substantial fraction of the AGB wind, and blows two opposite jets that shape the 
nebula into a bipolar structure. Point symmetric structures are also possible due to precession 
of the accretion disk around the companion. 

(3) Elliptical PN with large departure from sphericity are formed by stellar companions in 
a common envelope. The companion in the envelope does not blow strong jets (or blows 
no jets at all), but the interaction with the envelope ensures high equatorial mass loss that 
leads to an expanding ring structure. 

(4) Elliptical PN with small departure from sphericity are formed by a substellar companion 
(a brown dwarf or a planet) that spins-up the envelope of the AGB progenitor. Weak jets 
may be present. 

The framework above was based primarily on a series of key observational results and 
in part on a handful of theoretical studies, which we summarise below: 



The theoretical arguments was then, and is still now, that the axi-symmetric structure 
(including; po i nt-symmetric ) requires the presence of a stellar or substellar companion 
(ISokerl l2001bl . 120041 . l2006bl ). The rarer point-symmetric PN cannot be explained by a 
model based on shaping by the fast wind during the PN phase, and seem to require 
precessing accretion disks. Theoretical models showed th at binary inte ractions have 



the ability to produce the observed PN morphologies (e.g.. iMorrisI 119871 ) . 



• From population synthesis studies the fraction of PN that were derived from binary 
interaction was far lower than the needed 90% (that is the entire PN population minus 
10% spherical PN). For instance, lHan et al.l (119951 ) derived that only 38 ± 4%, of all 
PN derive from binary interactions (common envelopes, mergers and wider binary 
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interactions). 



We knew observationally that ~ 10% of all PN are round. We also knew that ~ 15% 
of PN are bipolar while the remaining ~75% are elliptical (mildly, extremely, and with 
or without addi tional structures). The d etailed classification does change from author 
to author (e.g. 
subdivisions. 



Corradi fe Schwarzlll995[ ) so we refrain here from describing additional 



At the time it appeared that the PN obser yed around p ost-common envelope central 
stars were not bipolar, but rather elliptical ( ISokerlll997| ). 



Bipolar PN have relatively higher expansion velo cities, more in line with the escape 
velocities of main sequence rather than AGB stars ( ICorradi fe Schwarzlll995l ). Elliptical 
PN, on the the hand, have expansion velocities that can be explained by winds from 
AGB stars. 

The spherical halos of many PN, are much fainter than the elliptical inner region. In 
the binary model this is accounted for by a binary interaction that occurs at a late 
stage in the superwind phase, the perio d of enhanced mass- loss that is observed to 
take place towards the end of the AGB (IDelfosse et al.lll997l ). The interaction both 
increases the mass- loss rate an d causes the departure from axisymmetry; even a planet 
can play a role in this process (jSokerl 120001 ) . 



Planet statistics were poor in that year. 



Later, Soker (2001a), realising the possible importance of wider stellar companions, 
added a fifth class of wide (but not very wide) companions that accrete mass and blow jets. 
The shaping is only by the jets. A moderately elliptical PN is formed, with pronounced 
signat ure o f jets (although not a s strong as in scenario 2 above). Finally, |Pe Marco fe Moe 
( 120051 ) and lSoker fe Subagl ( 120051 ) introduced the concept that not all 1 — 8 M stars actually 
make a PN. Single star s can m ake a spherical PN of a sub-luminous nature that is harder to 
detect; ISoker fe Subagl ( 120051 ) termed this class "hidden PN". 

In Table [TJ we summarise the conclusions of Soker (1997; column 3), Soker (2001a,b; 
column 4) and Soker & Subag (2005; column 5), based on the results above, as well as 
the conclusions of the present work, that will be justified in the following sections. Once 
the spherical PN are accounted for by mechanisms involving no interaction and bipolar PN 
are accounted for by close companions outside the common envelope (~10% and ~15%, 
respec t ively) , the remaining ~75% of all PN needs to be accounted for by other interactions. 



Sokerl (119971 ) expected that the only other PN shaping interaction would be a common 
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envelope interacti on, and thes e interactions are not frequent enough to account for 75% of 
all PN. This lead ISokerl ( 119971 ) to conclude that there are not enough stellar companions to 
produce the large fraction of non-spheric al PN. Even including much wider binaries in the 
class of interacting binaries (ISokerll2001a( ). there was still a need for a relatively large number 
of planetary companions to help shape PN. 



3. In the light of new results 

Since the mid-2000s, several results both of a theoretical and observational nature, 
allowed us to refine the conclusions of Soker (1997,2001a,b) and Soker & Subag (2005). Here 
we divide the new results into two classes. Those pertaining PN and those that regard the 
presence and frequency of planets around main sequence and evolved stars. 



3.1. Highlights of relevant PN research in the last decade 



• The fact th at shaping axi-symmetric PN with single stars may be difficult was discussed 
before (e.g.. lMorrislll987i : ISoker fc Liviolll989l: ISoker!ll997i l2006al: iNordhaus et alil2~0~0~7 



but see also De Marco (2009)). Recently (ISokerl l2006bl ; INordhaus et~al1l2007h it has 
become even more implausible. Global magnetic fields in AGB stars tend to feed back 
negatively onto the star's spin, slow it down and shut down the magnetic field itself in 
less than ~100 years, i.e., shorter than it is needed to sha pe the superwind >1 00 yr . 
Magnetic fields in the successful shaping model s of. e.g.. iGarcia-Segura et al. fiooj 
are imposed and constant. The suggestion of iBlackman et al.l (120011 ) that tapping 
part of the convection energy to sustain the magnetic -field forming dynamo is a viable 
mechanism to sustain AGB fields, w as questioned by ISoker &: Zoabil (120021 ) as well as, 
later on, by the authors themselves (INordhaus et al.l 120071 ). 



• Several theoretical papers have studied the effects of a binary companion on a mass- 
losing AGB star. These papers show that the relation between the binary interaction 
type and resulting PN mo r pholo gy may not be as str aight forward as envisaged by 
the simple scheme of ISokerl (119971 ): the seminal work of iMastrodemos fc Morris! (119991 ) 
studied the effects of a comp anion accreting from an AG B st ar wind and forming a 
stable accretion disk and jets, iReyes-Ruiz fc Lopea (119991 ) and INordhaus fc Blackman 
(120061 ) studied the formation of accretion disks and jet s around the primary's core by 
disrupted companion. Further, IBlackman et al.l (120011 ) concluded that more massive 
companions may blow powerful jets that can innate bipolar lobes. Equatorial flows can 
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arise from these systems (jSoker 



be formed by c ompanions in comm on envelope interactions (jSandquist et al.lll998l ) and 
in wider orbits ( jEdgar et al.ll2008r). bu t it is not clear whether bipolar morphologies can 



19971 ) . Finally the PN morphol ogy may change after it 



is ini tially established at the hand of the fast post-AGB wind (IHuarte Espinosa et al. 



20101 ) . although it is not known exactly which types of structures are susceptible to 
this change and which are not. 

Population synthe sis simulations t hat predict the fraction of PN shaped by binary 



interactions (e.g., lHan et al.l Il995l ). start from the assumption that all stars in the 
main sequence m ass range 1 — 8 will f orm a PN. Using a different population syn- 



thesis technique, iMoe fc De Marcol (120061 ) predicted the number of PN (with radius 
< 0.9 pc) present today in the Galaxy, if all of the 1 — 8 M star s actually make a 
PN (4 6 000=bl3 000 objects). This number was recently updated by lMoe &: De Marco 
( I2010l ) to 61 000±17000 objects. They compared this prediction with the actual num- 
ber of observed PN (with radius <0.9 pc - this radius limit is key to avoid much larger 
uncertainties that would derive from PN detection) of 8000±2000, derived from extra- 



galactic PN counts ( ]Jacobylll980l ). or 11000±2000, derived from a local sample count 
( iFrewl 120081 ). The predicted number is discrepant with the observationally-based es- 
timate of the Galac tic PN population at the 3a level. This calculation supports the 
earlier suggestion of ISoker fc Subagl ( 120051 ) that not all stars in the 1 — 8 M range ac- 
tually make a PN. If so, the PN population does not need to mirror the main sequence 
population, but it derives from only a subset of it. In Table [1] we reserved a row to 
"naked central stars". We exchange the term hidden PN used by Soker & Subag (2005) 
with naked central star of PN, to reflect the observational consequence that these will 
not be classified as a PN at all. Still, we do predict that very deep observation may 
reveal very intrinsically faint PN that are expected to be spherical. While Soker & 
Subag (2005) listed together the spherical PN, whether observed or hidden, in Table 1 
we classified them into two groups: spherical PN and naked central stars. 



Ziilstral ( l200lluDe Marcol (120091 ) and lMiszalski et al.l (12009bl ) catalogued the morpholo- 
gies of PN associated with post-common envelope binaries, concluding that there is a 
preferential association between these binaries and bipolar PN (although not all post- 
common envelope central stars do have a bipolar PN). This releases the constraint 
of Soker (1997) to associate bipolar PN with binary systems that avoided a common 
envelope interaction. In recent years the issue of what conditions dict ate that a binary 
avoids a common envelope interaction h a s been hotly debated (e.g., iBeer et al.l 120071 ; 



Carlberg et al.ll200a Bear k Sokerll2010l ; iNordhaus et al. 



2010), including when a bi- 



nary enters a common envelope at a late stage ( ICorradi et al.ll201ll ). Therefore, at this 
point we consider all close binary interactions (common envelope and not) in the same 
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group. 



• The MASH survey JParker et all l2006bl fl: [ Miszalski et al.ll2008l ) revised the fraction 
of spherical PN from ~10% to ~19%, where the fraction of spherical PN is higher 
for fainter/older PN. Decreasing the fraction of non-spherical PN to ~81%, this result 
releases slightly the pressure on the needed number of binary interactions. The MASH 
survey also slightly decreased the fraction of bipolar PN to 13%0. 

• Further, if very wide stellar companions shape mass-loss by accreting AGB wind and 
blowing jets, the number of possible binary interactions that can create non spher- 
ical PN is greatly increased. An example of such wide interactions is Mira AB, 
where despit e the large orbital se paration of ~ 40 times the AGB stellar radius 
(338-402 R ; IWoodruff et al.l |2004| ), the companion influences the wind morphology 



Matthews fc Karovska 2006; Karovska 



in the region between the two stars (IMarengo et al.ll200ll ; iKarovska et al.lll997l . 12005 



200J). 



Miszalski et al.l (j2009al ) revised the fraction of PN binaries with perio d shor t er tha n 



~15 days to 17±5% (only slightly larger than the 10-15% determined by lBondl (120001 )). 
This fraction is a lower limit, due to the photometric precision of the survey used, 
and the various survey biases (such as the angle between the orbital axis and the 
line of sight). How e ver, a ccounting for such biases would not increase this fraction 
dramatically: iBondl (120001 ) had assumed that the period distribution of the detected 
systems (P<3 days) meant that longer periods could not be observed and so a large 
fraction of binary PN may be hidin g at periods slightly longer than 3 days. On the 



other hand iMiszalski et al.l (j2009al ) determined observationally and |Pe Marco et al. 



(120081 ) determined theoretically that the detection limit of this survey technique is ~2 
weeks and that the paucity of central star binaries in the period range 3 days to 2 
weeks is real. We can therefore state that while the fraction of post common envelope 
binaries is only slightly larger than 17%, the fraction of central stars in binaries with 
periods longer than 2 weeks remains unknown and could be large. 

• We still do not know the fraction of PN which have binary central stars in the period 
range 2 weeks to the max imum period that still allows for an interaction. Preliminary 
work by iDe Marcol (120101 . see also Passy et al. (2011, in preparation)) indicates that 



1 We note here that of all the MASH PN 14% are either stellar (unresolved), irregular or asymmetric. 
The remaining PN are ellipticals (54%), round (19%) or bipolar (13%). In this paper we absorb the 14% of 
stellar, irregular or asymmetric PN in the elliptical class, such that elliptical, round and bipolar PN comprise 
68%, 19% and 13%, respectively. Clearly this is likely to artificially inflate the elliptical class, but since this 
is the most numerous class the percentage error brought by this inclusion remains relatively small. 
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this fraction may be larger than predicted by the current PN scenario where almost all 
l-to-8 stars make a P N. If confirmed, this findin g is in line with the conclusions of 



Soker & Subad (120051 ) and Moe fc De Marcol (l2006h . 



3.2. Highlights of relevant exoplanet research in the last decade 



About a decade ago estimates of the percentage of planet-hosting solar-like stars stood 
at ~ 3%, although it was already suspected that the percentage might be much higher 
for stars with hi gher metallicities: ~ 25% — 30% for stars w ith twice the solar metallicity 
(Fe/H > 0.3; e.g. lsantos e7aDl2004j : [Fischer fc Valentill2005h . iLineweaver fc Gretheri (120031 ) 
extrapolated from the detected parameter space of planetary systems to below detection 
sensitivity and concluded that the real fraction of planet hosting Sun-like stars should be at 
least 9% for M p smi > 0.3Mj and P < 13 years and at least 22% for M p sini > O.lMj and 
P < 60 yrs. They also suggested that since this area of the planet mass-period plane is only 
a fraction of that occupied by our own Solar System, these fraction may be even larger once 
the entire parameter space is sampled. This was the situation when Soker & Subag (2005) 
updated the PN formation channel statistics (See Column 5 of Table 1). 

Today we have strong evid ence that the planet ary fraction increases with both metal- 



licity and mass of the host star (IJohnson et al. 



20101 ) . although our knowledge of the planet- 



hosting star fraction as a function of pla net mass and orbital separation is grossly incomplete. 
A new finding in this respect is that of iBowler et al.l ( 20101 ) who determined that ~ 26^g% 



of stars having a main sequence mass of 1.5 < M/M < 2.0 host massive planets at large 
separations (but still less than 3 AU). The new finding is not only of a larger fraction of 
planet-hosting stars, but also puts the planets at larger orbital separations, where they can 
more easily be available to shape winds from AGB stars: if the planets are too close to their 
host, they interact with the star during the RGB phase, and either prevent the star from 
reaching t he AGB at all (and form a PN), or they just get destr oyed before the star reaches 
the AGB iNelemans fc Taurisl[l998l : Ivillaver & Liviol[20071 . l2009h . There is also an indication 
from the analysis of a handful of thick disk stars that the planetary fraction for low metal 



licity, old stars may be higher than for younger stars at similar metallicities (ISheehan et al. 
2010h . 



Another relevant finding is that of brown dwarf and planetary mass companions around 
extreme horizontal branch (EHB) stars, also termed sdO and sdB subdwarfs. Although such 
findings regard stars that have gone through the RGB, not the AGB evolution, we can draw 
information of a general interest. EHB stars are on the horizontal branch ( HB) and have 



small envelope mass, because their RGB progenitor lost most of its envelope (ID 'Cruz et al. 
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19961 ). Current consensus is that in most c ases the mass los s was caused or enhanced by the 
interaction with a stellar- mass companion (IHan et al.ll2002l ). 



few AU from single ( 


Silvotti et al.ll2007) 


2011; 


Qian et al. 


2009a 


bj) EHB stars. 1 



3) EHB stars. For the latter group, the companions orbit close 



binaries that went thought a common envelope interaction when the RGB progenitor 
of the sdB star engulfed its stellar companion. It is not clear where these tertiary/low 
mass companions were located within the pre-common envelope binary, but they may 
have survived despite the dynamical mayhem of the common envelope interaction^. 



Geier et al.l (120091 ) discovered a companion with mass 8 — 23Mj (either a planet or 
a low mass brown dwarf) with period of only 2.391 days around HD149382. This 
is the first detection of a substellar object, with a lower possible mass reaching the 
planet domain, that definitely went through a common envelope evolution inside an 
RGB envelope. The survival of the planet implies that interactions between gas giants 
and stars alter stellar evolution, therefore if a planet at a suitable distance is present 
around a growing AGB star, a common envelope with such planet may be survived 
and lead to shaping o f the ejecta. A second possible planet (Msini = 1.25 Mj) was 
recently discovered by lSetiawan et al.l (120101 ) at 0.116 AU from the metal poor HB star 
HIP 13044. 



4. Updating the PN shaping framework 



In light of the ne w discoveries 



U) a 



w e reframe here the res ults of ISokerl (119971 ) and 



subsequent updates by ISokerl (l2001al ibl) and ISoker fc Subaa (120051 ). The work discussed in 



§ 13. II releases the tight constraint that brought ISokerl (119971 ) to conclude that a large number 
of PN should be shaped by planets. The increasing numbers of brown dowarf and planet 
discoveries around evolved stars also give support to the idea that planets play a role in 
shaping PN, provided that they are at a suitable distance from the star to interact during 
the AGB. Finally, the larger fraction of planets detected recently around main sequence 
stars, compared to what was believed a decade ago, provides us with more flexibility when 
discussing the role of planets in shaping PN. 

In what follows we quote percentage figures with an accuracy of 1%. This accuracy 
is unrealistically high. However, we do so for ease of following some of the arguments. At 



It is not excluded that these planets formed in the ejected common envelope (|Peretsll2010l) 
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the end, we will discuss what should be considered reasonable errors on these percentage 
estimates. 

We start by returning to the discussion of § 13.11 instead of considering the PN population 
as the inevitable child of the 1 — 8 M star population, needing to reflect its binary fraction, 
we can think of the PN population as deriving from only a subset of the 1 — 8 M stars. 



We adopt the subdivisions of PN morphological types envisaged by ISokerl (119971 ). with the 
changes brought in by the MASH morphological discoveries (§|5J Table [2]). 

First of all, we determine the f raction of all 1 — 8 M stars that result in a PN of 
any shape. iMoe fc De Marcol feoioh concluded that only 18±8% of all 1 - 8 M Q stars 
(11000/61000) make PN and 82±8% make no PN, or do not go through the AGB at all 
(see § I3.ip . This estimate rests on single stellar evolution, the initial mass function and 
the galactic star formation history, as well as observations. It is almost independent of the 
binary fraction and period distribution. 

An independent way to determine the fraction of all 1 — 8 M stars that result in a PN 
of any shape is to consider only the main sequence binary fraction and period distribution. A 
summary of how these fractions are determined can be found in Table Approximately 57% 
of F and G main sequence st ars are in binaries and ~30% have an orbital separation < 30 AU 
( IDuquennoy &: Mayor! Il99ll ). and will interact sometime during the primary star's life. We 
have taken a maximum separation that is possibly too large for some types of interaction, 
but is smaller than the distance between Mira and its interacting companion (see § 13. ip . 
We assume, once again, that an interaction is needed to produce non-spherical PN. We also 
account for the fact that 1 4% of all stellar systems su ffer a strong interaction on the RGB 
(orbital separation < 3 AU; IDuquennoy fc Mayorlll99ll ) that induces so much mass-loss that 
the system is prevented from ever evolving to the AGB. We then deduce that 30 — 14=16% 
of all stars suffer an AGB interaction with a stellar companion and go on to form a non 
spherical PN. This means that ~16% of all stars go through an AGB interaction, ~14% go 
through and RGB interaction and never ascend the AGB and the remaining ~70% suffer no 
interaction with a stellar companion and evolve into a naked central stars or a central star 
with round PN, unless they have an interaction with a substellar companion. As lSokerl (119971 ) 
we assume that mildly elliptical PN are shaped by interactions with substellar companions. 
Hence, ~16% of all stars go through an interaction with a stellar companion and result in 
non-round and non-mildly elliptical PN. This is equivalent to stating that ~16% of all stars 
result in 61% of all PN (the non-round and non mildly elliptical ones; Table [3]), and this 
means that 26% (16/61x100) of all 1 — 8 M stars make a PN of any shape. This estimate, 
which has a probable error of ±5%, is based on approximate binary con siderations, and is 
consistent with that outlined above (cf. 26% with 18±8%), obtained by IMoe &: De Marco 
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( 120 lOl ) almost entirely from single stellar evolution arguments and PN counts. 



In Table S] we finally adopt 26% as the fraction of 1 — 8 M stars that are able to make 
a PN of any shape. Approximately 14% is the fraction of all stars that do not ascend the 
AGB because of a strong interaction on the RGB. The remaining ~60% is the fraction of all 
stars that do not make a PN because they suffer no interaction (for further discussion see 
§ [5]). Finally, we split 16% into fractions reflecting the PN morphological types presumed 
to derive from binary interactions with stellar companions, i.e., 13% are bipolar (16x0.13 / 
0.61=3.4), 28% are extremely elliptical (16x0.28 / 0.61=7.3), and 20% are mildly elliptical 
with jets (16x0.20 / 0.61 =5.2). The remaining 10% of all 1-8 M Q stars (26-16%) goes to 
form the two remaining PN classes, the 19% round PN (10x0.19 / 0.39=4.9) and the 20% 
mildly elliptical PN with no jets (10x0.20/0.39=5.1). 



Recently iMoe &: De Marcol (120 lOf ) carried out a more precise binary population syn- 
thesis model, completely independent of morphological subdivisions, and determined that 
8% of all 1 — 8 M Q stars suffer a strong interaction with a stellar companion on the AGB 
(common envelopes and interactions outside the envelope). This is in line with the approxi- 
mate estimate based on binary considerations and morphological subtypes, where ~11% of 
all 1 — 8 M stars suffer a strong AGB interaction (3.4+7.3%: Table H] and last column of 
rows 2 and 3 in Table [T|) . 

In Table [T] (last column) we finally list a revised subdivision of the evolutionary channels 
followed by intermediate mass stars and their link to PN morphology. During the last 
decade of simulations and observations we have learned that a one-to-one correspondence 
between morphologies and shaping channels is as improbable as the likelihood of finding 
two identical PN. This is why we stress (see also note "c" in Table [1]) that these are only 
guidelines. The take-home messages are that (i) common envelope and other interactions, 
where the companion is close to the AGB surface, likely result in strong collimation leading 
to bipolarity and extreme ellipticity; (ii) that low mass companions such as planets are likely 
only to promote small departures from a spherical shape; (iii) that blowing jets may be 
promoted not only by an accreting stellar companion at some distance from the primary, 
but also by a disrupted lower mass companion that forms a disk around the primary core 
( INordhaus fc Blackmanll2006l ). If so, this would confuse the distinction between channels 4 
and 5 in Table [TJ 



If, as in ISokerl ( 119971 . l2001al ) and ISoker fc Subagl ( 120051 ) , we assume that mildly elliptical 
PN with no jets were shaped mainly by a substellar companion, we predict that substellar 
companions have shaped ~20% of all PN, which corresponds to ~5% of all 1 — 8-M stars 
having suffered an interaction with a planetary companion on the AGB. However, if substellar 
companions can get disrupted, form a disk around the primary core and blow jets, then this 



percentages could rise to as much as ~40% and ~8%, respectively. 



4.1. The implied fraction of planets from the AGB perspective 



The fact that only ~30% of main sequence stars have stellar companions close enough 
to interact, combined with the fact that ~81% of PN have non-spherical shapes that appear 



to need a companion, lead ISokerl ( 119961 ) to state that "substellar objects (brown dwarfs or 
gas-giant planets) are commonly present within several AU around main-sequence stars. For 
a substellar object to have a high probability of being present within this orbital radius, 
on average several substellar objects must be present around most main-sequence stars of 
masses < 5M ." 

We have derived here that ~5% of all 1 — 8 M stars suffer a strong interaction with 
a substellar companion on the AGB. This estimate depends critically on (i) the assumption 
that mildly elliptical PN with no jets are mostly shaped by a substellar companions, (ii) 
on the frequency (~20%) of this specific PN morphology and (iii) on the two independent 
arguments that predict that that only 18% or 26% of 1 — 8 M stars make a PN at all (0.18 
or 0.26x0.20 ~ 0.05). 

These substellar companions should occupy orbits between ~3 and ~30 AU around the 
main sequence progenitors. Closer companions (whether stellar or substellar) suffer a strong 
interaction (likely a common envelope interaction) on the RGB and are either destroyed or 



preclude the primary from ascending the AGB because of excessive mass- loss ( iSoker et al. 



1998 ; Bear fc Soker 2010 ; Nordhaus et al. 2010 ). Companions farther on% would not inter- 
act. 

Using a population synthesis technique one could use this estimate and work backwards 
to determine the parent population of these star-planet systems. Even without such cal- 
culation, we can already draw a few conclusions. The median main sequence progenitor 



mass of today's PIN 



(IMoe & De Marco 



popu lation is 1.2 M and the median metallicity is approximately solar 



2Q06J). Takin g the maximu m radii on the AGB, Ra, and on the RGB, 



Rr, from llben fc Tutukovl (119851 ). ISokerl ( 119981 ) derived the following approximation for the 



3 The limit out to which a co mpanion is captured tidally during AGB phase, according to Soker ( 1996h . is 
^12 AU. iNordhaus et al.l (|2010l ) would revise this limit down. We have however widened it to 30 AU in light 
of other types of interactions that do not involve captures. We have cited the AGB star Mira as an example 
of a wide binary which is nevertheless interacting. Whether this limit applies to low mass companions and 
planets is not known. 
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Table 1. Evolutionary channels for different morphological types. 





Evolutionary channel 


By-product 


1997 


Percentage 
2001 


2005 


This work 


1 


No interaction 


Spherical PN 


~ 10 


~ 10 


~ 10 


~5 a (19) b 


2 


Close stellar companion 
outside envelope 12 


Bipolar PN 


~ 11 


~ 15 


~ 15 


~ 3 a (13) 6 


3 


Stellar companion in 
a common envelope 


Extremely 
elliptical PN 


~ 23 


~ 25 


~ 25 


~ 7 a (28) 6 


4 


Substellar companion 
in a common envelope'' 


Mildly 
elliptical PN 


~ 56 


~ 35 


~ 15 


~ 5 a (20) 6 


5 


Wide stellar 
companion d 


Elliptical PN 
+jets 




~ 15 


~ 15 


~ 5 a (20) 6 


6 


No interaction 


Naked central 
stars (post-AGB) 






~20 


~ 60 e 


7 


Strong interaction 
on the RGB 


EHB stars 
(post-RGB) 








~ 14 e 



Note. — Subdivision of the evolutionary channels taken by all 1 — 8 Mq stars and the evolution 
of our understanding o f the percentages over the past decade. The percentage are from Soker 
(19971. ISokeil l|2001al lbh. lSoker fc Subael {200d \ and the present paper. EHB stars are formed form 
RGB stars that have lost most of their envelope due to a strong binary interaction with a stellar 
companion and will never ascend the AGB. 

a All these percentages should be considered "a few percent" . However, their relative accuracy 
may be much better than their absolute one. They add up to 25 instead of 26 because of rounding 
errors. 

b All numbers in brackets are fractions of the PN population and add up to 100. The distribution 
of percentages within brackets reflect the distribution of morphological types. 

c We maintain these two channels separate for historical reasons, although we know that these 
two evolutionary channels do not always result in the corresponding morphologies. They could be 
combined into one channel called: strong binary interactions. 

d It is conceivable that also substellar companions may blow jets if they ar e destroyed in a common 
evnelope interaction and form a disk around the primary's core (§ 13. H and [Nord haus fc Blackmail 

l2poeh . 

c These uncertainties are of the order of ±5%. 
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Table 2. Percentages of all PN that reside in different morphological classes. 



Percentage of PN 



Class 



Reference 



81% 
19% 
13% 
28% 
20% 
20% 



non-round PN 
round PN 
bipolar PN 

extremely elliptical PN 

mildly elliptical PN with no jets 

mildly elliptical PN with jets 



Parker et al. 2006b 

Parker et al. 2006b 

Parker et al. 2006b 

Parker et al. 2006b : Soker 1997° 

Parker et al. 2006b : Soker 1997° 

Parker et al. 2006b: Soker 1997° 



a Sixty-eight percent elliptical flParker et al.l l2006bl) = 28% + 20% + 20% The subdi- 
yision into three classes follows ISokerl (119971 ). The jet /no jet division is from lBalick et al 
dl998h . 



Table 3. Subdivision of all 1-8 M Q stars into orbital separation groups. 



Percentage 


Class 


Method of 


of 1-8 M stars 




determination 


57% 


binaries 


DM91 a 


43% 


single stars 


100 - 57% 


30% 


binaries with a < 30 AU - interaction on the RGB or AGB 


DM91 a 


70% 


single stars and binaries with a >30 AU - no interaction 


100 - 30% 


14% 


binaries with a < 3 AU - interaction on the RGB 


DM91 a 


16% 


binaries with 3 < a < 30 AU - interaction on the AGB 


30 - 14% 



Duquennoy &: Mayor! 1 1991 
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Table 4. Subdivision of all 1-8 M stars into binarity and interaction classes. 



Percentage 


Class 


Method of 


of all 1-8 M stars 




determination 


(14 a +e b )% 


do not ascend the AGB and do not make a PN 


DM91; Table E 


3% 


make bipolar PN 


16% c x0.13/0.61 d 


7% 


make extremely elliptical PN 


16% c x0.28/0.61 d 


5% 


make mildly elliptical PN with jets 


16% c x0.20/0.61 d 


5% 


make spherical PN 


(26-16)% c xO. 1970.39 


5% 


make mildly elliptical PN 


(26-16)% c x0.2070.39 


26% 


make PN of any shape 


~ (3 + 7 + 5 + 5 + 5)% 


(60-e b )% 


ascend the AGB but do not make a PN 


(100 - 26 - 14)% 



a This fraction is tied to the choice of the orbital separation limit (we have selected 3 AU; see 
Table [3]) within which binaries and star-planet systems interact on the RGB. This limit is not 
well known. It also depends sensitively on the primary mass and the strength of the tides. 

b e denotes a small fraction of stars that have an interaction with a substellar companion on 
the RGB and, as a result, does not ascend the AGB. The same small fraction is taken out of the 
group of stars that do ascend the AGB. 

c Sixteen percent is the percentage of all stars that have a stellar companion and that suffer an 
interaction on the AGB, resulting in bipolar, extremely elliptical and mildly elliptical with jets 
PN; 26% is the percentage of all stars that make a PN of any shape, so 26-16 is the percentage 
of all stars that make round and mildly elliptical PN; see Table |3j 

d These are the morphological percentages from Table [2J 
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ratio of the maximum AGB to maximum RGB radius as a function of mass: 

\og(R A /R R ) = 3.71og 2 (M/M ) - O.371og(M/M ) + 0.16, M < 2.25M , (1) 

for stars which develop degenerate helium cores, and 

\og(R A /R R ) = 2.2 - 1.81og(M/M ), M > 2.35M , (2) 

for more massive stars, where M is the primary's mass on the zero age main sequence. From 
these relations we can state that for stars with spectral type earlier than A the radius ratio 
is > 2, while for stars with spectral type later than F-G, which correspond to the progenitors 
of central stars of PN, this ratio is < 2. For the progenitors of PN, therefore, the planetary 
population has to be relatively far out. We have quoted limits of 3 and 30 AU, respectively. 
he lower limit depends critically on the adopted prescription o f tidal capture and mass-loss 



flSokerlll996l : IVillaver fc Livid 120071 . 120091 ; iNordhaus et al.ll2010f ) 



Considering that the fraction of solar-like stars hosting relatively close by planets may 
be of the order of 10-20% and considering that we predict that a low ~5% of all solar like 
stars have substellar companions further out than ~3 AU, it is probable that the substellar 
companions farther out are the outermost companions in a system of multiple planets. 



5. Do we observe the naked central stars? 



In this subsection, we want to elaborate further on the suggestion that ~60% of all 
1 — 8 M stars that ascend the AGB do not make a visible PN. This is a surprisingly high 
number and implies a large population of post- AGB stars which, under the commonly as- 
sumed scenario, would have P N, but in the revised scenario do not. To be more precise, 
Moe fc De Marcol (120061 . 120101 ) predict that ~61000 stars in the Galaxy are in the evolu- 
tionary phase appropriate to be surrounded by a PN of less than 0.9 pc in radius, i.e., these 
stars are in a post- AGB with T eS > 30 000 K, have left the AGB less than 25 000 years ago, 
and have mass > 0.55 M . This population is about 6 times larger than the Galactic PN 
population with radius smaller than 0.9 pc (11000 PN; iFrewl 120081 ) . from which we predict 
that 50 000 (61 000 - 11 000) stars in the Galaxy are "naked" central stars. It would be much 
harder to detect naked central stars, because of high reddening in the Galactic plane. How- 
ever, we should con sider that the deep MASH survey has now detecte d ~3000 PN, which is 
<30% of the total jjacobv! fl98ol : iMoe fc De Marcol I2OO6I : iFrewlhoOsh . PN are particularly 
bright and have emission lines in the red part of the spectrum which are less subject to 
reddening. The fraction of detectable naked central stars should be far smaller. 



Naked central stars would look like subdwarf O and B stars, most of which are post- 
RGB stars, and would be easily confused with them. In some cases they could be told apart 
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from those sdOB stars which are in a post-RGB phase of evolution on the grounds of higher 
luminosity, and larger surface gravity (alt hough there is an overlap between post-RGB and 
post-AGB stars on the log g — T e ff plane (lNapiwotzkilll999l their Fig. 4)). 



The fraction of naked central stars in the sdOB sample should be only a few percent 
of the total, because of their short lifetimes compared to the post-RGB phase. Any post- 
AGB object, independent of its mass, will fade below 100 L in less than 100,000 years 
JVassiliadis & Woodfl994T ). All post-AGB objects with M > 0.60 M will fade to L = 100 L 
i n only 10,000 years. HB lifetimes are instead of the order of several tens of millions of years 
( IDorman et al.lll993f ). 



We have an indication that the predicted number of post-AGB sdOB stars are there, but 
this jnust be corroborat ed by a more accurate count. O'Tool (private communication - but see 
also iHirsch et all 120081 ). determined that in about 120 sdO stars, there are several objects 
which could be in the post-AGB phase. This could constitute the few percent predicted. 
For the sdB stars, a few hundred of which have determine d logq and T e ff, s everal reside 
in th e logg — T e ff locus appropriate for post-AGB stars ( iNapiwotzkil Il999t IHirsch et al. 



20081 ). There should be fewer sdB post-AGB stars than sdO post-AGB stars, because their 
evolutionary times for the temperature T e ff < 40 000 K is quite a lot shorter than for 
T eff >40 000 K. 



Looking at Fig. 2 of lNapiwotzkil ( 119991 ) or Fig. 3 of IHirsch et all ( 120081 ) we see that there 
are a number of sdB and sdO stars that are not surrounded by a PN. Several objects are 
indeed found on the horizontal part of the post-AGB tracks, where a PN is expected. Others 
are seen on the descending part, but at high core mass such that one might expect that a 
PN would be present, considering the shorter evolutionary timescales. This sample is not 
homogeneous, and it is therefore difficult to determine whet her the number of naked centra l 
stars on this plot is what is expected from the prediction by lMoe fc De Marcol (120061 . 120101 ) . 



Deep imaging surveyes may reveal a number of faint and spherical PN around naked 
central stars. The fact that the deep MASH survey approximately doubled the fraction of 
circular PN may already have bourn out this prediction. 



6. Summary 

The first assumption we adopt in this paper, which is increasingly supported by obser- 
vations and theoretical considerations, is that single stars are mostly unable to produce PN 
whose shapes diverge from spherical. Second, we bring to bear the population prediction 
that the fraction of stars in the initial mass range 1 — 8 M that actually form PN is only 
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^20% (see § 3.1). 



As in the past we make the simplistic assumption that mildly elliptical PN with no 
jets are formed exclusively by planetary interactions. This assumption is broadly justified 
by (i) simulations that show that density contrasts drastically reduce for dec reasing com- 
panion mass in commo n envelopes (ISandauist et al.lll998l ; iDe Marco et al.ll2003l ) and outside 
the common envelope ( Kim fc Taam 201ol ); such interactions are therefore likely to lead to 
mildly elliptical shapes at best, (ii) The fact that known PN around close binaries are by and 
large bipolar or strongly elliptical and have often jets and substructures; if the stellar-mass 
companions are at larger orbital separations than those that interact in a common envelope, 
theoretical considerations suggest that accretion onto the companion promotes jets, which, 
once again, lead to more dramatic departures from spherical symmetry than a mild ellip- 
ticity. Despite these justifications, the one-to-one correspondence between interactions with 
substellar companion and the generation of mildly elliptical PN should be considered only 
as a guideline to be confronted by observations, and against which to continue testing future 
shaping theories. 

With the premise above, we predict that the fraction of PN shaped by planets (~ 20%) 
corresponds to only ~ 5% of all 1 — 8 M stars having interacted with a planet on the AGB, 
indicating that a few percent of 1 — 8 M stars should have Jupiter-class companions further 
out than a few AU, a thing that se ems increasingly plausible given new disc overies of planets 



around main sequence stars (e.g., iJohnson et al.ll2010l ; iBowler et al.ll2010l . see section 3.2 



Planets farther out may be the outer planets in planetary systems. 

Finally, the implication that ~ 60% of all stars in the initial mass range 1 — 8 M Q do 
not go through a PN phase, may be justified observationally. There are several hot subdwarf 
stars that occupy a location of the log g — T eS diagram appropriate for central stars of PN, 
but that do not exhibit a PN. In order to strengthen this claim one would have to make a 
prediction of the fraction of all subdwarf O and B stars that are expected to be post-AGB 
in origin through a population synthesis. Also, one would have to homogeneously locate a 
sufficient number of these stars on the log g — T e g diagram to determine if the prediction is 
verified. 
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